Disproportionation Transition at Critical Interaction Strength: Nai/ 2 Co02 
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Charge disproportionation (CD) and spin differentiation in Na!/ 2 Co02 are studied using the 
correlated band theory approach. The simultaneous CD and gap opening seen previously is fol- 
lowed through a first order charge disproportionation transition 2Co 3 ' 5+ — > Co 3+ +Co 4+ , whose 
ionic identities are connected more closely to spin (S=0, S=| respectively) than to real charge. 
Disproportionation in the Co a g orbital is compensated by opposing charge rearrangement in other 
3d orbitals. At the transition large and opposing discontinuities in the (all-electron) kinetic and 
potential energies are slightly more than balanced by a gain in correlation energy. The CD state 
is compared to characteristics of the observed charge-ordered insulating phase in Na!/ 2 Co02, sug- 
gesting the Coulomb repulsion value U is concentration-dependent, with U(x — 1/2) ~3.5 eV. 

PACS numbers: 71.28. +d,71.27.+a,75. 25. +z 



The discovery by Takada and coworkers of 
superconductivity^] in Na .3CoO 2 »1.3 H 2 near 5 
K has led to extensive studies of the rich variation 
of properties in the Na^CoOaiSystem (0.2 < x < 1). 
It has been pointed out@, M that the formal 
charge of the Co ions that occupy a triangular lat- 
tice, whose simple average is 4 — x, changes at 
x ~ 0.5: for x > 0.5 the observed Curie-Weiss (C- 
W) susceptibility indicates (1 — x) Co 4+ ions with 
spin S=^ and x Co 3+ with S=0, while for x < 0.5 
the susceptibility is Pauli-like and all Co ions are 
identical ("Co( 4 - x )+"). In a narrow region around 
x = 0.5 there is charge ordering (CO) as observed 
by electron diffraction, 0, U |(J probably accompa- 
nied by antiferromagnetic spin ordering. While the 
rest of the phase diagram is metallic, this x = 0.5 
phase undergoes a CO+metal-insulator transition 
(MIT) at 50 K as verified by resistivity^ and optical 
conductivity. 0, 01 

The mechanism of charge ordering has been a cen- 
tral question in many classes of transition metal ox- 
ides, especially the cuprates and manganates. The 
CO question is beginning to be addressed in these 
cobaltates, from the mechanism of ordering at com- 
mensurate concentrations 0, 0, to the effects of 
fluctuations when the concentration varies from a 
commensurate fraction. [2j, llCj Preceding the ques- 
tion of CO however must come the more basic one of 
disproportionation of a lattice of identical Co' 3+a: ' + 
ions for x < 0.5 to disproportionated Co 3+ /Co 4+ 
ions with their distinct charges and spins. CD is 
a signature of intraatomic correlation and is closely 
associated with local moment formation, but it also 
has a collective nature to it (every ion dispropor- 
tionates). Moment formation has been studied pri- 
marily in single band models in conjunction with 
the correlation-driven MIT, [llj but CD in multiband 
systems (as Na .5CoO 2 almost certainly is) is a sub- 



stantially more involved question and is relatively 
unstudied. 

Recently some of the present authors obtained a 
CD/MIT0 in a study of Na x Co0 2 (x = 1/3, also 
x = 2/3) using the correlated band theory method 
LDA+U (local density approximation plus Hubbard 
U). For no apparent reason, the CD/CO transition 
appeared simultaneously with the MIT (band gap 
opening) as the interaction strength U was varied in 
the 3-4 eV range. Another feature that was not clari- 
fied is that this CD involved rather little actual redis- 
tribution of charge but at x = | showed up dramati- 
cally in the magnetic moment as the metamorphosis 
3Co 3 f + -► Co 3+ +2Co 4+ . One of the central ques- 
tions in these cobaltates is the relevance of its multi- 
band nature, a complicating aspect that has stimu- 
lated studies into the possibility of orbital-selective 
metal-insulator transitions. [Tj, Ei ^ n Na a: Co0 2 the 
threefold i 2ff manifold is partially filled, with bro- 
~~* a g + e 'g arising from the Co0 2 
LDA calculations indicate both 
i g auu o_ j.-Gi.jjm aurfaces already at x = |, and the 
presence of e' g Fermi surfaces has become a central 
component of several models of the superconducting 
phase, [il) 

Another important question is that of corre- 
lation effects and especially the mechanism of 
CD/CO/MIT transitions. There have been a few 
pioneering applications of correlated band theories 
to model disproportionation. 0, 0, ^1 In this paper 
we study more thoroughly the x — | case, which 
has become more compelling since a CO/MIT was 
observed. Q The variable in our study will be U, but 
in a small parameter range such as near the critical 
value U c , varying U will be analogous to applying 
pressure to change the U/W ratio, hence predictions 
are open to direct experimental comparison. In ad- 
dition, the CD/CO that is observed ass^ 0.5 sug- 



ken symmetry i 2g 
layered structure. 

and e' Fermi surfaces already at x 




FIG. 1: Effect of U on the magnetic moments for 
Na^CoCh. The jump reflects a first-order transition (and 
accompanying hysteresis) in the critical region (3.2 eV to 
3.6 eV). Gap opening (metal-insulator transition) simul- 
taneously occurs with CD due to the first-order nature. 
Dashed and solid lines represent the CD and the undis- 
proportionated state, respectively. Insert: Change of 
the magnetic moments for the O ions: Ol (with 2n site 
symmetry) is shared by two Col and one Co2, while 02 
(with 2m site symmetry) is shared by one Col and two 
Co2. 



gest a concentration dependent interaction strength 
U(x). We study in detail a mean-field but fully self- 
consistent treatment as provided by LDA+U, of a 
first order CD transition driven by correlation, and 
show that the multiband nature of Na^CoC^ is an 
integral ingredient of the CD transition. We com- 
pare with data to suggest the value of U(x). 

Our calculations that are based on the underly- 
ing hexagonal structure (space group P6$22) hav- 
ing lattice constants (a = 2.84 A, c — 10.81 A).[j"t| 
were done with a two Co supercell with space group 
P2/m) allowing two ions Col and Co2 to be real- 
ized in the cell. The O height was taken as zq = 
0.168(c/2) = 0.908 A as relaxed by LDA calcula- 
tion by Singh. [13] The calculations were carried out 
with the full-potential nonorthogonal local-orbital 
minimum-basis scheme (FPLO) [l8| and both pop- 
ular schemes [l9ll20T | for LDA+U functional with the 
main results being common to each. Since Na a; Co02 
is a good metal (except at x — \ which we address) , 
the results quoted in this paper are those from the 
so-called "around mean field" scheme 19] appropri- 
ate for relatively small U. The Brillouin zone was 
sampled with 128 (306 for LDA+U calculations) k 
points in the irreducible wedge. Orbitals for the ba- 
sis set contained 2s2p3s3p3d for Co, O, and Na pro- 
vide excellent basis flexibility. 

The doubled unit cell allows for CD into Co 3+ and 
Co 4+ , but with only one magnetic ion there is no is- 
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FIG. 2: Change with U of d charge for the two Co ions. 
Upper panel: change of total 3d charge (the charge is 
transferred to O ions). Lower panel: changes of the ma- 
jority and minority separately. 



sue of spin order (except simple alignment). Hence 
we work throughout with a ferromagnetic (FM) dou- 
bled cell. In the local spin density approximation a 
FM state is favored (by a small energy) over the 
nonmagnetic state, with a total moment of 1 /is as 
noted first by Singh. ^3] Because the Na ions are 
ordered above the CO transition, [2l| the Na site is 
nonsymmetric with respect to the Co sites, and the 
Co ions have somewhat different moments because 
the symmetry is already broken by the Na ordering. 

The effect of the on-site Coulomb repulsion U on 
the magnetic moments is displayed in Fig. ^ We 
follow behavior with U continuously by using the 
solution from one value of U as the starting point 
for the case with incrementally increased/decreased 
U. Although unconstrained, the total cell moment 
remains at 1 fig: half metallic in the metallic phase, 
insulating in the CD phase. Both Co magnetic mo- 
ments decrease slowly (-0.025 /is/eV) with increas- 
ing U, probably due to some charge transfer Co — > 
O. At U=U C 2=3.6 eV, the solution changes character 
discontinuously, with the small and large moments 
providing the clear identification Col— >Co 3+ , Co2-^ 
Co 4+ . By U — 4.5 eV the large moment obtains a 
maximum and the small moment is negligible. For 
larger U the Co 4+ moment tends to decrease. Upon 
decreasing U the moments vary continuously down 
to U = U c i = 3.2 eV, below which the solution jumps 
back to the initial undisproportionated state. Hence 
we find a 0.4 eV hysteresis in this CD transition 
centered at U c =3.4 eV. The first order nature is 
consistent with no symmetry having been broken 
in our choice of simulation cell, and also accounts 



2 



llF 



0.8 



5 0.6 



a 
s 



0.4 



0.2 - 



iiYiYiYt =6= = --^r~"^--^ 




□ — □ Col (up) 
A— A Co2 (up) 
0—0 Col (dn) 
V— 7 Co2 (dn) 



4 5 

U(eV) 




FIG. 3: Change of the occupancy n ag of the a 9 state 
versus U, which reveals the strong a g charge dispropor- 
tionation at the transition. Majority orbitals of both 
Co ions are fully occupied regardless of U . At U = U c , 
the occupancy of the minority decreases for Co2 and in- 
creases for Col by 0.23 e. 



FIG. 4: Change in energy E(U) per doubled supercell 
versus U, relative to E(U=3 eV), showing the various 
contributions discussed in the text. Values shown are 
energy differences with respect to the value at U = 3 
eV, except the LDA+U energy Ea whose actual value is 
plotted. Insert: Blowup of the critical region. 



for the simultaneous CD and gap opening observed 
previously jj|: the regime of accelerating dispropor- 
tionation includes the opening of the gap, but this 
region is inaccessible because of its higher energy 
than either the undisproportionated or CD states. 
This first order transition is analogous to low-spin 
«-> high-spin transitions mapped out using fixed spin 
moment calculations. |22| but generalizing to a two 
dimensional space of Col and Co2 moments 
respectively as considered first by Moruzzi, Marcus 
and Kiibler.[23 

Although Col and Co2 become quite different 
above U c , the disproportionation does not show up 
as strongly in Co 4+ — > Co 3+ charge transfer as in 
the spin rearrangement. As U increases from zero, 
a small amount of charge from both Co ions trans- 
fers (0.04-0.05 e/eV) onto the O ions (this shift is 
smaller if using the "fully localized limit" LDA+U 
functional 20]). At the transition, as shown in Fig. 
121 the Col^Co 3+ charge is continuous at the tran- 
sition, while the Co2^Co 4+ charge jumps by ^0.1 
electron. The discontinuities in the majority and 
minority charges separately are also pictured in Fig. 

m 

The picture of the charge rearrangement at the 
transition becomes clearer when the a g occupation 
is monitored, as shown in Fig. [3] At the CD tran- 
sition, much of the minority a g charge on Col (— > 
Co 3+ ) transfers to the majority on Co2 (becoming 
Co 4+ ). The total charge difference is only a frac- 
tion of the change in minority a g charge, because 



it is compensated rather strongly by a rehybridiza- 
tion and back polarization of the other 3d {e' g and 
e g ) states, reminiscent of what was seen for varia- 
tion of the Na content. [24| A large change in spin 
is achieved without large change in total charge be- 
cause of the multiband nature of the system at and 
near the Fermi level. 

The energy E(U) is shown in Fig. 21 I n the 
critical region, the undisproportionated state is fa- 
vored marginally (only 8 meV/Co), a difference so 
small that the total energy is essentially continuous 
at the transition. As CD occurs, the (Kohn-Sham) 
kinetic energy jumps sharply by 0.74 eV/Co while 
other contributions decrease: (in eV/Co) 0.46 for 
the (electron+nuclear) potential energy, 0.09 for the 
LDA exchange, and 0.18 for the LDA+U correlation 
energy (cx U). The LDA+U correction is modest but 
essential to the transition, and only becomes large 
enough to tip the balance at U — U c . 

Now we discuss the relation of our study to 
Na^OoC^. Na orders in this system above the 
MITQ so our Na-ordered cell is realistic in that re- 
spect. The observed cell is twice as large as we have 
allowed, so it is not possible to study AFM order- 
ing in our cell. In x = | samples, Co 3+ and Co 4+ 
ions are disproportionated (C-W susceptibility) but 
charge/spin disordered above the transition; conduc- 
tion is poor and abruptly becomes insulating when 
CO occurs. In our (T=0) simulation, Col and Co2 
are undisproportionated with a half metallic spec- 
trum for U < U c , whereas for U > U c the ions dis- 
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proportionate. Within LDA+U, the MIT is orbital- 
selective, pja [T3'J While the distinction between our 
results and cobaltate behavior are evident, the ques- 
tion of CD is quite separate from that of CO: experi- 
mentally, the C-W susceptibility indicates CD in the 
sample at high temperature, but the simultaneous 
CO/MIT occurs only at 50 K. 

From the infrared conductivity Hwang et al. have 
reportedQ a very small 10 meV gap in their in- 
sulating x — i sample. This small value suggests 
the value of U required to model the x = ^ sys- 
tem is very close to (only slightly larger than) U c , 
i.e. U(x = i) ~3.5-4 eV. From susceptibility, heat 
capacity, and resistivity data the £=0.7-0.8 phases 
are clearly more correlated (U > 4 eV), confirmed 
by the very large field-dependent thermopower [25) . 
an extremely large Kadowaki- Woods ratio ,2(|, and 
many reports of magnetic ordering. On the other 
hand, in the superconducting phase x ~ 0.3 the re- 
sistivity is lower and the mass enhancement over the 
band value is very small 2, 4] so this phase is at most 
weakly correlated with U(x ~ 0.3) < 3 eV. 

Whereas previous LDA+U studiesJ^J of the 
Na x Co02 system did not obtain disproportionation, 



we have shown that reducing constraints allows, 
and leads to, disproportionation. In opposition to 
LDA+U conventional wisdom ( "precise value of U is 
not very important"), results for Na x Co02 are very 
sensitive to the value of U around the critical value 
U c = 3.4 eV. We view NazCo02 in terms of a con- 
centration dependent U(x), which may be in part 
a guise for a single band — * three band crossover 
around discussed earlier. Q In such cases 

studies 28] of the multiband Hubbard model suggest 
U e ff for the single band regime (x close to unity) 
becomes U e ff = U c j\fi for the three band regime 
(x < i), in which case U for the superconducting 
phase could be even smaller than our estimate above. 
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